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ABSTRACT 

Motivated by recent discoveries of low-density super-Earths with short orbital periods, we 
have investigated in-situ accretion of H-He atmospheres on rocky bodies embedded in dissipating 
warm disks, by simulating quasi-static evolution of atmospheres that connect to the ambient disk. 
We have found that the atmospheric evolution has two distinctly different outcomes, depending 
on the rocky body's mass: While the atmospheres on massive rocky bodies undergo runaway disk- 
gas accretion, those on light rocky bodies undergo significant erosion during disk dispersal. In 
the atmospheric erosion, the heat content of the rocky body that was previously neglected plays 
an important role. We have also realized that the atmospheric mass is rather sensitive to disk 
temperature in the mass range of interest in this study. Our theory is applied to recently-detected 
super-Earths orbiting Kepler- 11 to examine the possibility that the planets are rock-dominated 
ones with relatively thick H-He atmospheres. The application suggests that the in-situ formation 
of the relatively thick H-He atmospheres inferred by structure modeling is possible only under 
restricted conditions; namely, relatively slow disk dissipation and/or cool environments. This 
study demonstrates that low-density super-Earths provide important clues to understanding of 
planetary accretion and disk evolution. 

Subject headings: planets and satellites: formation 



1. INTRODUCTION 

An increasing number of low-mass exoplanets 
with masses of less than about 20 M^, which 
arc often called super- Earths (SEsfl, have been 
detected recently, thanks to progresses in radial- 
velocimetry and transit photometry including 
operation of two space telescopes, Kepler and 
CoRoT. The size measurement allows us to access 



^ While exoplanets with masses of 10-20 are sometimes 
called exo-Noptunes instead of super-Earths, we simply call 
them super-Earths in this paper. 



the planetary interiors theoretically and thereby 
estimate the planetary bulk compositions. It has 
been revealed that there are many low-density SEs 
that are larger in size than rocky (i.e., iron-silicate) 
objects of the same masses. As for SEs, the vari- 
ety of possible ingredients obscures the composi- 



composition (Valencia et al. 


2007: Seasjer et al. 


2007: Adams ct al. 


2008: 


Grasset et al.l 20091 


Roeers & Seaecri 


201C 


). In particular, low-density 



SEs are subject to the degeneracy. Sometimes, the 
degeneracy has a crucial impact on understanding 
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the origin of planets and planetary systems. 

A typical example is degeneracy originating 
from the uncertainty of the presence and mass 
of H-He atmospheres. The mass-radius relation- 
ship for a planet alone is insufficient to distin- 
guish whether the planet is a rock-dominated or 
water-dominated one, because the H-He atmo- 
sphere makes up the difference in size between 
rocky and water bodies of the same mass. Since 
transiting SEs in general are orbiting close to their 
host stars (namely, well inside the snow line), 
whether they are rocky or water planets would af- 
fect considerably understanding of the accretion 
and migration history of the planets. An example 
of that is recently-detected SEs orbiting a Sun-like 
star named Kepler-11. 

The Kepler-11 system is a multiple-planet sys- 
tem that contains a t least five low-density SEs 
(jLissauer et al.ll201l[ ). Of them, the densities of 



Kepler-lld and lie are as low as 0.9 ^[j g gem ^ 
and 0.5 gcm~^, respectively, which are lower 
than those of pure-water planets. This suggests 
the presence of H-He atmospheres on the planets 



Indeed , the structure modeling bv iLissauer et al 



poill ) suggests a possibility that those two SEs 
are rock-dominated planets with thick H-He atmo- 
spheres that account for 10-20 % of the planetary 
masses. If so, they are a new type of planet that 
we have never seen in the current solar system. 

The H-He atmospheres, if primordial, came 
from the protoplanetary disk where the planets 
formed. In the framework of the core accretion 



model (jHavashi et al.l 119851) . solid planetary em- 



bryos first form, and then collect the ambient disk 
gas gravitationally to form H-He atmospheres. If 
this proceeds well before disk dispersal, the disk- 
gas accretion enters a runaway state at some crit- 
ical point, which results in forming envelopes of 
gas g iants (|Mizunolll980t iBodenheimer &: PoUackl 
19861 ) which are much more massive than the pre- 
dicted atmospheres of the Kepler-11 planets. The 
critical point is when the atmospheric mass is 
appro ximately 1/4-1/3 of the protopl anet's total 
mass (Stevenson 1982 : Wuchterl 19931) — the frac- 
tion being similar to those of the atmospheres of 
interest. This means that the inferred masses of 
the H-He atmospheres of the Kepler-11 planets are 
close to critical, which motivates us to investigate 
the possible amounts of H-He gas that proto-SEs 
with short orbital periods gain from protoplane- 



tary disks. 

In this study, we explore the possibility of the 
in-situ accretion of H-He atmospheres. As for the 
origin of H-He atmospheres of short-period SEs, 
possibilities of Neptune-like planets (i.e., water- 
dominated SEs) th at form far out an d migrate 
close to host stars (jRogers et al.ll2011l ) and rem- 
nants of gas giants whose envelopes were stripped 
(jNavakshin .201 L) were previously discussed. The 
in-situ accretion of H-He atmospheres on rocky 
SEs from disks has not been explored by direct 
simulation of disk gas accretion in the mass and 
orbital period regimes of transiting SEs. The in- 
situ accumulation of H-He atmospheres would be 
reasonable in the modern picture of planet for- 
mation. Short-period SEs are likely to have mi- 
grated to their current locations. The promising 
mechanism to move them inward is the type-I mi- 
gration winch occurs via planet-disk tidal inter- 
action ( Ward Il986h . The protoplanets formed in 
this way are packed closely together due to reso- 
nance c apture, according to recent iV -body simu- 
lations ( Terquem &: Papaloizoul 2007 ). After the 
migration, the disk dissipates. In many cases, 
the decline of disk gas density triggers dynami- 
cal instability of a multiple-protoplanet system, 
resulting in orbital crossing and collisio ns of the 
protoplanets (e.g., Ogihara fc Idal 20091 ). Subse- 
quently, the protoplanets collect the surrounding 
disk gas to form atmospheres in the dissipating 
disk. Thus, the competition between the atmo- 
spheric accumulation and the disk dissipation may 
yield intermediate-mass atmospheres like the pre- 
dicted atmospheres for the Kepler-11 planets. 

This paper is organized as follows. In sec- 
tion [21 we describe our theoretical model of the 
atmospheric accumulation with emphasis on ef- 
fects that we newly incorporate in this study. 
Then, we show numerical results of the atmo- 
spheric growth and the sensitivities of the final 
atmospheric masses to parameters involved in the 
model. In section |3l the masses of the accreted 
atmospheres that we calculate are compared with 
those of the atmospheres of Kepler-11 planets in- 
ferred by internal-structure modeling. Based on 
the comparison, we discuss the possibility of the 
in-situ accretion of the atmospheres. Other pos- 
sibilities are also discussed. Finally, we conclude 
this paper in section |4l 
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2. ACCUMULATION OF ATMOSPHERE 
2.1. Model 

Wc simulate the radially ID structure and 
quasi-static evolution of the atmosphere of a pro- 
toplanet embedded in a protoplanetary disk. The 
atmosphere's contraction (or expansion) results in 
its mass gain ( or loss). The de t ail of the model 
is described in llkoma fc Genda The nu- 

merical integration is done with the code that w e 
developed and used in iHori fc Ikoi^(|20ld[201ll ). 
Below is a brief summary of the model. 

The planet consists of a compressive atmo- 
sphere with solar elementary abundance on top 
of an incompressive "rocky" body with density 
of 3.9 g/cm^. Choice of value of the rocky den- 
sity and incorporation of rocky small compress- 
ibility have tiny impacts on the structure and 
mass of the atmosphere. We integrate a usual 
set of four equations describing the spherically- 
symmetric, quasi-hydrostatic structure of a self- 
gravitating atmosphere, which includes the equa- 
tions of hydrostatic equilibrium, mass conserva- 
tion, radiative/convective energy transfer, and en- 
ergy conservation (i.e., time change in entropy). 
The equation of state for th e atmospheric gas tha t 
we use in this study is from lSaumon et alT ( 19951 ). 

The atmosphere is assumed to be equilibrated 
with the disk at the smaller of the Bondi and Hill 



spheres fsee llkoma fc Gendal 120061 ): namely, the 
atmospheric density and temperature are equal to 
those of the ambient disk gas there. On the other 
hand, the bottom of the atmosphere corresponds 
to the interface between the atmosphere and the 
rocky body. The atmosphere is heated by the un- 
derlying rocky body at the bottom. In this study, 
the energy flux from the rocky body is given at the 
atmospheric bottom as an inner boundary condi- 
tion, as described below. 

In this paper, we consider grain-free atmo- 
spheres, unless otherwise noted, to investigate 
upper limits to the masses of the H-He atmo- 
spheres that the SEs gain via the in-situ accretion. 
Thus, we assume that the atmospheric opacity in- 
cludes only the gas op acity, which is taken from 



Freedman et al.l ()2008l) in this study. We suppose 



that the protoplanet formed via giant impacts in 
a dissipating disk. Each simulation of the atmo- 
spheric growth starts with an arbitrarily hot state 



(i.e., a given high luminosity). The initial choices 
of the luminosity do not affect our results because 
the initial cooling of the atmosphere occurs fast 
relative to the disk dissipation. No continuous en- 
ergy supply by planetesimal bombardment is as- 
sumed. This assumption also leads to finding up- 
per limits of the atmospheric mass. 

In this study, wc include two effects newly: We 
assume that disk dissipation occurs concurrently 
with the atmospheric growth. Furthermore, we 
incorporate the effect of the heat capacity of the 
rocky body. While both effects were not included 
in previous studies, they have crucial impacts on 
atmospheric growth in the situations considered in 
this study, as shown below. 

The disk gas density, pd{a, t), is assumed to de- 
crease in an exponential fashion, namely. 



Pd{a, t) = p"{a) cxp {-t/Td), 



(1) 



where a is the semi-major axis, t is time, and 
is the dissipation time that is regarded as a pa- 
rame ter: the i nitial disk density, taken 
from iHavashil ( 1981 ) in this study. According to 



modern theories of disk evolution (Calvet et al 



2000l for a review), the overall dissipation first 
occurs via viscous diffusion in ~ lO'^ yr. Once the 
gas density becomes so low in the inner disk that 
the stellar UV penetrates through the disk, photo- 
evaporation occurs around the gravitational radius 
of several AU. From that time on, the inner disk 
is separated from the outer disk and evolves via 
viscous diffusion separately. Since this is equiva- 
lent to the fact that the disk size decreases by a 
factor of 10, the diffusion timescale (being propor- 
tional to the square of disk size) for the inner disk 
decreases by a factor of 100 to become ^ 10^ yr. 
Based on a simple analytical argument, namely, 
solving the ID viscous-diffusion equation for sur- 
face gas density, S^, 



dt 



a da 



(2) 



with = at the inner and outer edges, one 
finds that the inner disk dissipates in such a way 
that pd (X exp(— TT^i/lO^yr), which corresponds to 
Td ^ lO^yr. This is consistent with results from re- 
cent numerical si mulations of disk evolution (e.g., 
Gorti et aLlliooj ). 

Decrease in di sk density results in co oling of 
the atmosphere ( Ikoma fc Genda 20061 ). The 
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rocky body below the atmosphere also cools 
down. Detailed treatment of the rocky body's 
thermal evolution is beyond the scope of this pa- 
per. Instead, we evaluate its contribution by in- 
putting the luminosity from the rocky body in 

the form of Lj-ock = - Crock A^rock^Tl/dt + -^radio 

as an inner boundary condition for the atmo- 
spheric structure, where Crock is the specific heat 
of silicate (= 1.2 x 10^ ergR-^g-^), Mrock the 
mass of the rocky body, Tt the temperature at 
the atmospheric bottom, and Lradio the luminos- 
ity due to the radioactive decay of chondrites 



2 X 10^°(Mrock/Mm)ergs~^ (see iGuillot et al 



19951 ). In the present case, since the decline of 
Tfc is caused by that of pd, 



dTb _ dlnpd dTb 



1 dn 



dt d In Pd Td d In pd 



(3) 



Furthermore, from the analytical solution o f 
the radiative atmosphere ( Ikoma fc Genda 20061) . 
Tf, is related to the disk temperature, T^, as 
dTb/dlapd — Td/i. Hence, we calculate Lrock 



as 



ock 



-^^rockCrockTd 
4Td 



-^radi( 



(4) 



2.2. Atmospheric growth 



Examples of atmospheric growth are shown in 
Figure [1] In the simulations, Mrock is assumed 
to be 4 M^. The blue line shows the atmo- 
spheric growth when the disk gas density is con- 
stant through the simulation. In this case, after 
gradual growth, runaway gas accretion starts at 
t ~ 3 Myr, resulting in forming an atmosphere 
that is much more massive than the atmospheres 
of interest in this study. 

The green and red lines show the atmospheric 
mass evolution with concurrent disk dissipation 
with Td = 1 X 10^ years; the former and latter 
represent the cases with Crock = (i.e., no contri- 
bution of the rocky body's cooling to heating the 
atmosphere) and Crock = 1-2 x lO'' ergK~^g~^, 
respectively. As shown by the two lines, the at- 
mospheric growth levels off at some point in these 
cases, because the disk is assumed to disperse be- 
fore the runaway gas accretion would happen if 
it were not for disk dispersal. After the leveling- 
off, the atmosphere is found to be eroded. The 
simulations were stopped when the disk gas den- 
sity became low enough that the mean-free path 




X 



IC 10" 

Time [yr] 

Fig. 1. — Evolution of the atmosphere on the 4Mg3 rocky 
body for disk temperature = 550 K. In the blue-hne case, 
the disk density, p^j, is assumed to be constant through the 
simulation. In the green-line and red-line cases, the disk 
is assumed to dissipate as oc exp (— t/lO^yr), which is 
shown by a thin gray line. Also, the rocky heat capacity, 
Crocki is 1.2 X 10^ ergK~^g~^ in the red-line case, while 
Crock = in the green-line case. Note that the radioactive 
luminosity is also included, but has a negligible impact on 
the atmospheric evolution. 



of the ambient gas molecules was longer than the 
planetary radius; 



kTd 



Pd 



1.4 X 10 



-19 



Mr 



ock 



IMe 

Mrock 



Td \ , 3 /r\ 

g/cm (5) 



100 K 



IM, 



ViooKy V 



lAU 



11/4 



where k is the Boltzmann constant, G the gravita- 
tional constant, and a the coUisional cross-section 
of the molecule (= 2.5 x 10"^^ cm^). 

The atmospheric erosion during disk dissipa- 
tion has been newly found in this study. The 
erosion is due to atmospheric expansion. The 
expansion occurs, because disk depressurization 
enhances the outward pressure gradient near the 
outer boundary, which pushes the atmosphere out- 
wards. Then, the atmospheric gas seeps out from 
the protoplanet's gravitational sphere, which con- 
tinues until the atmosphere is equilibrated with 
the ambient, deprcssurized disk. The decrease in 
atmospheric mass due to this effect (the green line) 
is, however, small, relative to the decrease due to 
heatup by the cooling rocky body (the red line). 
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< 0.001 




Mass of rocky body [Earth mass] 

Fig. 2. — The final atmospheric mass as a function of the 
rocky body's mass. The red solid line represents a fiducial 
model where the disk temperature, Tj, is 550 K, the disk 
dissipation time, r^j, is 100 kyr, and the initial disk density, 
p^, is equal to that of the minimum-mass solar nebula, 
Pmsn llHavashilll98ll ): the values of the three parameters 
are used below unless otherwise specified. The red dashed 
and dotted lines are for t^j = 10 kyr and 1 Myr, respec- 
tively. The green and blue lines are for T^j = 940 K and 
200 K, respectively. In the ycUow-line case, is 10 times 
as large as Pmsn- While grain-free opacity is assumed in 
the other cases, the grain opacity for protoplanetary disks 
from lSemenov et al. ] ||2003| ') is used in the purple-line case. 



As described in section 12.11 the disk dissipation 
cools the surface of the rocky body, resulting in 
heat supply to the atmosphere. This lifts up the 
atmosphere, which results in further atmospheric 
erosion. This is why the atmospheric erosion is 
more significant when the rocky heat capacity is 
incorporated. 

2.3. Final masses of the H-He atmo- 
spheres 

In Fig. [21 the final mass of the H-He atmo- 
sphere, Mn+He, which is divided by the rocky 
body's mass, Afrock, is shown as a function of 
Afrock for different values of the disk's tempera- 
ture, Tdi and dissipation time, Td- While those 
two parameters are focused on in this study be- 
cause upper limits of Mn-i-He are of special inter- 
est, the sensitivity of A'/h+Hc to the initial disk 
density, p^, and the atmospheric opacity is briefly 
seen below. 

First, for a given set of the parameters, Mh-i-ho 
increases with Mrock in the small-Mpock regime. 
Also, it is found that there is a critical value of the 



rocky body's mass, M*^^^, beyond which the atmo- 
sphere becomes quite massive eventually. This is 
because the atmospheric accretion enters the run- 
away state before the disk is depleted significantly. 
Thus, the atmospheric growth has two distinctly 
different outcomes, depending on A/rock- The at- 
mospheric growth beyond the critical point is dis- 
cussed in section [3l 

The three red lines of different types show the 
sensitivity of A/h+hc to r^. For a given Afrock, 
larger results in larger Afjj-i-Hc- The sensitivity 
is not large; Afn+Ho increases by a factor of less 
than 10 for two-order-of- magnitude increase in t^. 
Also, M*Qj,]^ increases, as decreases. This is be- 
cause a short lifetime of the disk (i.e., small Td) 
allows only massive protoplanets to start runaway 
accretion. Difference in Td also makes only a small 
change in M*^^^^ and also A/h-i-Hc at the critical 
point (< O.lA/iock)- It is noted that while physi- 
cally small, the sensitivities of A/jj-f He may suffice 
to validate observationally inferred masses of the 
atmospheres, as discussed in section [3l 

The results for different three disk tempera- 
tures are shown by the solid green {Td = 940 K), 
red {Td = 550 K), and blue {Td = 200 K) with 
Td = 100 kyr in Fig. H It is found that Td 
has a significant impact on Afn+Hc and M*^^^^. 
Qualitatively, high Td results in small A/h-i-Hc 
and, therefore, large M*^^^. This is basically be- 
cause atmospheres made from warmer gas (more 
exactly, higher-entropy gas) are gravitationally 
less bound. However, the impact of the outer 
boundary conditions on the atmospheric mass is 
known to be negligibly small in the case of massive 
atmospheres embedded in cool disks like proto- 
envelopes of gas giants. This is because most of 
the atmospheric mass is concentrated in the deep 
region of the atmosphere whose structure is insen- 



sitive to the outer boundary conditions ( Mizuno 



19801 : LStevenscmlll982t IWuchterllll99,l llkoma et al 



20011 ) ■ In contrast, in the case of low-mass atmo- 



spheres embedded in warm disks which are ap- 
propriate to SEs with short orbital periods, our 
results in Fig. [2] reveal that the impact is signifi- 
cantly large. 

The initial disk density, p'^, has only a small im- 
pact on A/h+Hc and M*^^^^. For example, shown 
by the yellow line is the result for the case of p° 
being 10 times larger than that used for the red 
solid line. Since high density corresponds to low 
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entropy, the atmosphere embedded in denser disks 
tends to be more massive for the same reason de- 
scribed above. However, it is realized that the at- 
mosphere before being eroded is massive enough 
that the structure in its deep part is insensitive to 
disk properties. 

Finally, the values of Mh+hc and M*^^^^ that 
we have derived are upper and lower limits, re- 
spectively, for a given set of and in the sense 
that no grain opacity and no additional energy in- 
put are assumed. In the case represented by the 
purple line, the effect of grain o pacity is evalu- 
ated. To do so, we have adopted ISemenov et al 



(|2003[ )'s protoplanetary-disk grain-opacity model. 
This would be an extreme case of large opacity, 
namely, a lower limit to A/h+Hc because large 
opacity sl ows down contraction of the atmosphere 
(jlkoma et al.,.20oA . Thus, the actual solutions are 
in between the lines of grain-free (red solid line) 
and grain-rich (purple line) atmospheres. 

3. APPLICATION TO KEPLER-11 SUPER- 
EARTHS 

3.1. Comparison with atmospheric masses 
inferred by modeling 

The super-Earths orbiting Kepler- 11 are thought 
to have relatively thick H-He atmospheres, as de- 
scribed in Introduction. In this section, we ap- 
ply our theory to the planets and, thereby, get 
some implications for their composition and ori- 
gin. Specifically, we compare the atmospheric 
masses derived from modeling with those which 
we have calculated in this study. 

The comparison is made in Fig. [31 The 
a;-axis is the planetary total mass, Mp (i.e., 
Mj-ock + .^H+He); the y-axis is the percentage of 
Mn+He relative to Mp. The black crosses repre- 
sent Mn+He derived by the structure modeling of 
the planets, which account for the observed masses 
and radii. The data were kindly provided by E. 
Lopez who calculated them using a modified ver- 
sion o f the modeling method from iLissauer et al.l 
(j201l[ ): the effect of the rocky heat capacity hav- 
ing been incorporated in the thermal evolution, 
which yielded a slight decrease in Mn+He (Lopez, 
personal communication). 

The curves show the results of our calculations 
in this study for four different cases. The temper- 
atures 940 K and 550 K correspond to those at 



100% 



10% 




0.1% 



Mrock+Mn+He [Earth mass] 



Fig. 3. — Application to the Kepler-11 planets. The 
black crosses represent the estimated atmospheric masses 
of Kepler-llb to llf, the data of which were provided by 
E. Lopez (personal communication). The red, green, and 
blue lines are for = 940 K, 550 K, and 200 K, respec- 
tively. The solid and dashed lines are for t^j = 1 Myr and 
10 kyr, respectively. We have drawn the extrapolated lines 
beyond the critical points represented by filled squares, as- 
suming that H-He gas accumulates on the rocky body with 
the critical mass. The points on the extrapolation lines are 
the results of simulations where the atmospheric accretion 
is limited by supply of gas from the disk that evolves via 
viscous diffusion; a- viscosity is adopted and the values of a 
are 1 X 10""' for = I Myr and 1 X 10"^ for = 10 kyr. 



Kepler-llb's and llf's semimajor axes (a = 0.088 
and 0.25AU) respe ctively, in an optically-thin 



disk (jHavashi Points indicated by filled 

squares are the critical points (see section 12. 3p . 
Beyond the critical points, we have also drawn 
extrapolated lines (hereafter called super-critical 
lines) by assuming that H-He atmospheres accrete 
on the non-growing rocky bodies with the criti- 
cal masses; mathematically, the super-critical lines 
are expressed by y = (1 — M*^^^/x) x 100%, where 
X = Afrock + Mh+ho and y = MH+Ho/(M.ock + 
-Mh+hc)- To check the validity of the super-critical 
lines, we have simulated the super-critical (i.e., 
runaway) accretion whose rate is not higher than 
supply limit of the disk gas due to viscous dif- 
fusion, SirTjdVd.-, where and Vd are the surface 
density and viscosity of the disk gas, respectively. 
We have adopted the a-prescription for the disk 
viscosity; = 10 kyr and 1 Myr correspond to 
a ~ 10~^ and lO"**, respectively. It is shown that 
the numerical solutions are certainly on the super- 
critical lines. In the case of = 550 K and = 
10 kyr, for example, Mrock = IS.OMg at the point 
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on the super-critical line, which is close to M*^^^. 
= 12.9 Mq; nevertheless, Mn+Ho/Mp is larger by 
a factor of ~5. This demonstrates that a slight 
difference in Mj-ock around M*^^^ results in a large 
difference in Mh+ho- 

Recent TV-body simulations of the f orma- 
tion of hot SEs (Terquem fc Papaloizoul 2007 : 
Qgihara fc Idal [looj" suggest that disk dissipa- 
tion often causes collisions between embryos of 
SEs. The major accretion of atmospheres occurs 
after the collisions. In the late stages of disk 
evolution, photo-evaporation separates the inner 
disk from the outer disk; then, the inner disk dis- 
sipates on a time scale of as short as 10^ years 



(|Gorti et al.ll2009l) . The dashed lines in Fig. H 



correspond to the cases of = 1 x 10 yr. Com- 
paring the lines with the crosses, one finds that the 
values of 7\/h+Hc inferred from modeling arc much 
larger than those derived in this study, except for 
Kepler-llb and 11c. Furthermore, because the at- 
mospheres are assumed to contain no grains, the 
derived atmospheric masses are maximal. Thus, 
the predicted atmospheres of the Kepler- 11 plan- 
ets are not accounted for by the in-situ accretion 
in such a quickly dissipating disk. 

Slower disk dissipation may be preferable in this 
respect. The green solid line shows the case of Td 
= 1 Myr and Td = 550 K. While the corresponding 
line for Td = 940 K is not drawn, the difference be- 
tween the solid and dashed lines is similar to that 
for Td ~ 550 K. In this case, as far as Kepler-llb to 
lie are concerned, the values of il/jj+Hc inferred 
from modeling seem to be consistent with those 
derived in this study. Kcplcr-llf lies above the 
line; only the high-mass end is close to the line. 
Furthermore, a cooler disk would be compatible 
with the Afn+Hc inferred for all the planets. As 
a reference, we show the result for Td = 200 K 
by a blue solid line. However, it would be worth 
noting that Kepler-lld and lie are near the crit- 
ical points in the super-critical regime. As seen 
above, such cases arc rare, because slightly large 
cores result in much more massive Affj+Hc- 

3.2. Mass Loss 

Before concluding this study, we discuss the loss 
of the atmospheres. The Kepler-11 planets are 
orbiting relatively close to their host star which 
is a G dwarf aged 8 ± 2 Gyr. Since the star is 
old, the current irradiation level of X-ray and EUV 
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Semimajor Axis [AU] 

Fig. 4. — The mass of atmospheric gas that a planet loses 
for 10 Gyr (solid lines). The numbers attached to lines in- 
dicate the estimated values. The mass loss is assumed to 
occur via the energy-limited hydrodynamic escape driven 
by XUV from a G dw arf like Kepl e r-11. Stellar-XUV data 
have been taken from lRibas et al ] l|2005l 1. Measured semi- 
major ax es and bulk densities of Kepler-llb, 11c, lid, lie, 
and llf l lLissauer et al.ll201lf) are shown by blue crosses. 
The numbers in parentheses are the likelihood values of 
the measured planetary masses. 



(XUV) is low. In the past, however, the planets 
should have been exposed to intense stellar XUV. 
Stellar-XUV irradiation results in loss of planetary 
atmospheres. 

Here we estimate the mass of the H-He atmo- 
sphere that a planet with a given density loses 
via the energy-limit ed hydrodynamic escape (e.g., 
Watson et al.lll98l[ ). We integrate the equation 



Ml. 



3efxu\ 



tide 



(6) 



where e is the heating efficiency, i^xuv the incident 
flux of stellar XUV, p the planetary bulk density, 
and i ^tirif- a correction factor for the Roche-lobe ef- 



fect (lErkaev et al.l2007tlLeca.velier des Etanes et al 



200. 



. We adopt e = 0.4 (see Valencia et al.ll201oir 



the fi tting formula of Fxuv from iRibas et al 
(|2005l ). and if tide = 0.92 appropriate for the 
Kepler-11 system. 

Figure |4] shows the mass that the planet loses 
for 10 Gyr as a function of semimajor axis and 
planetary bulk density, together with the mea- 
sured values of a and p for Kepler-llb to llf. 
It is revealed that despite of different distances 
to the host star, all the planets have lost simi- 
lar amounts of atmospheric gas, which are of the 
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order of 0.1 M®. Even for Kepler- lib {Mp = 
4.3+2-gMffi) and llf (Mp = 2.3±1.2Me) which are 
the two hghtest planets among them, the amounts 
of mass loss correspond to ~1-10 % of the plane- 
tary masses. (Although their early-stage inflated 
structure being neglected in this estimate, the re- 
sults are similar to those from more detailed sim- 
ulations by Lopez et al. [in preparation]). Thus, 
the consideration of mass loss does not affect the 
implications obtained above for Kepler-llc to llf. 
As for Kepler-llb, the condition becomes more se- 
vere. 

3.3. Other Possibilities 

3.3.1. Degassing 

Hydrogen-rich atmospheres may also be formed 
on rocky super-Earths via degassing. Plan- 
etesimals contain metallic iron, which is oxi- 



dized by water to produce hydrogen (jAbe et al 
However, detailed calculations show that 



the resultant hydrogen-rich atmosphere accounts 
for, a t most, several percents of the planetary 
mass ( Elkins-Tanton fc Seagerll2008t iRogers et al 



20111 ). This suggests that degassing is insufficient 
to explain, at least, the inferred masses of the 
atmospheres of Kepler-lld and lie. 

3.3.2. Water- dominated planets 

Another possibility is that those planets are 
water-dominated ones covered with relatively 
thin H-He atmospheres, namely, warm Neptunes, 
which would have formed beyond the snow line, 
followed by inward migration. The structure and 
formation of warm Neptunes were recently dis- 
cussed by [Rogers et al.l (j201l[ ). Looking at results 
of their modeling, one finds out that the atmo- 
sphere constituting < 1 % of the planetary mass 
is compatible with the observed radii of Kepler- 
lld and lie. Such atmospheres are stable for the 
current level of stellar XUV, as mentioned above 



(see also [Rogers et al.ll201l[ ). However, the origin 
is yet to be explored. Also, the resultant masses 
an d radii of the planet s formed in the simulations 
bv [Rogers et al. ( 2011 ) do not match those of the 
Kepler-11 planets, although they do not have a 
special focus on Kepler-11 in the paper. Since 
migration requires the persistence of the disk gas, 
further gas accretion should occur after arriving 
near the host star. Hence the origin of such atmo- 



spheres is also a challenging issue. 

4. CONCLUSION 

Motivated by recent discoveries of low-density 
super-Earths (SEs) with short orbital periods, 
we have investigated the in-situ accretion of H- 
He atmospheres on short-period SEs. Specifi- 
cally, we have simulated the quasi-static evolu- 
tion of the H-He atmospheres embedded in dis- 
sipating, warm disks and, thereby, derived the 
atmospheric masses at the time when disk gas 
has disappeared. We have also applied our the- 
ory to recently-detected short-period SEs orbiting 
Kcplcr-11, and examined the possibility that the 
planets are rocky ones with relatively thick H-He 
atmospheres. The application demonstrates that 
the in-situ formation of relatively thick H-He at- 
mospheres inferred by structure modeling is pos- 
sible only under restricted conditions; namely, rel- 
atively slow disk dissipation and/or cool environ- 
ments. 

Through the application of our theory to de- 
tected low-density SEs, we have gotten some im- 
portant suggestions for future improvement in the 
theory of the formation of H-He atmospheres of 
low-mass planets: (1) Atmospheric erosion occurs 
during disk dissipation. This suggests that not 
only when but also how protoplanetary disks dis- 
sipate are crucial factors for determining the at- 
mospheric mass. (2) The thermal contribution of 
the rocky body to the atmospheric evolution has 
a significant impact on the atmospheric erosion. 
Thus, it is important to investigate in more de- 
tail the thermal evolution of the rocky body that 
concurrently occurs with the accretion of the sur- 
rounding atmosphere. (3) In relatively hot en- 
vironments considered in this study, the atmo- 
spheric mass is rather sensitive to disk tempera- 
ture, which means that the late-stage thermal evo- 
lution of disks should be also taken into account 
in simulating the formation of H-He atmospheres 
of disk origin. 

Degeneracy in composition between rocky plan- 
ets with thick hydrogen-rich atmospheres and wa- 
ter planets is a critical issue in planet formation 
theories. Findings from this study are applicable 
to any short-period super-Earths and would be 
helpful for removing the degeneracy. As demon- 
strated in this paper, the atmospheric mass in- 
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ferred by structure modeling with observed mass 
and radius can be evaluated by comparison with 
that derived from the accretion theory. Also, low- 
density SEs provide important clues to under- 
standing of planetary accretion and disk evolution. 
Planetary candidates observed by Kepler contain 



abou t thousand SE-size objects (jBorucki et al 



I2OIII ). We expect that follow-up observations will 
identify them and determine their masses, which 
will lead to better understanding of the origins of 
low-mass planets. 

We are grateful to T. Guillot and J. Fortney 
for fruitful discussions. We also thank E. Lopez 
for kindly providing us with the data of his mod- 
eling of the structure and thermal evolution of the 
Kepler- 11 planets. We appreciate the anonymous 
referee's critical and constructive comments which 
helped us to improve this paper. M. I. is supported 
by Challenging Research Award from Tokyo Insti- 
tute of Technology and Grant-in- Aid for Scientific 
Research on Innovative Areas (No. 23103005) from 
the Ministry of Education, Culture, Sports, Sci- 
ence and Technology (MEXT) of Japan. Y. H. 
is supported by Grant-in-Aid for JSPS Fellows 
(No. 23003491) from MEXT, Japan. 

REFERENCES 

Abe, Y., Ohtani, E., Okuchi, T., Righter, K., & 
Drake, M. 2000, Origin of the Earth and Moon, 
413 

Adams. E. R., Seager, S., & Elkins-Tanton, L. 
2008, ApJ, 673, 1160 

Bodenheimer, P., & Pollack, J. B. 1986, Icarus, 
67, 391 

Borucki, W. J., Koch, D. G., Basri, G., ct al. 2011, 
ApJ, 736, 19 

Calvet, N., Hartmann, L., & Strom, S. E. 2000, 
Protostars and Planets IV, 377 

Cochran, W. D., Fabrycky, D. C, Torres, G., et 
al. 2011, ApJS, 197, 7 

Elkins-Tanton, L. T., & Seager, S. 2008, ApJ, 685, 
1237 

Erkaev, N. V., Kulikov, Y. N., Lammer, H., Sel- 
sis, F., Langmayr, D., Jaritz, G. F., & Biernat, 
H. K. 2007, A&A, 472, 329 



Frecdman, R. S., Marley, M. S., & Fodders, K. 

2008, ApJS, 174, 504 

Gorti, U., Dullemond, C. P., & HoUenbach, D. 

2009, ApJ, 705, 1237 

Grasset, O., Schneider, J., & Sotin, C. 2009, ApJ, 
693, 722 

Guillot, T., Chabrier, G., Gautier, D., & Morel, 
P. 1995, ApJ, 450, 463 

Hayashi, C. 1981, Progress of Theoretical Physics 
Supplement, 70, 35 

Hayashi, C, Nakazawa, K., & Nakagawa, Y. 1985, 
Protostars and Planets II, 1100 

Hori, Y., & Ikoma, M. 2010, ApJ, 714, 1343 

Hori, Y., & Ikoma, M. 2011, MNRAS, 416, 1419 

Ikoma, M., & Genda, H. 2006, ApJ, 648, 696 

Ikoma, M., Nakazawa, K., & Emori, H. 2000, ApJ, 
537, 1013 

Ikoma, M., Emori, H., & Nakazawa, K. 2001, ApJ, 
553, 999 

Lecavelier des Etangs, A., Vidal-Madjar, A., Mc- 
ConneU, J. C, & Hebrard, G. 2004, A&A, 418, 
LI 

Lissauer, J. J., et al. 2011, Nature, 470, 53 

Mizuno, H. 1980, Prog. Theor. Phys., 64, 544 

Nayakshin, S. 2011, MNRAS, 1274 

Ogihara, M., & Ida, S. 2009, ApJ, 699, 824 

Rogers, L. A., Bodenheimer, P., Lissauer, J. J., & 
Seager, S. 2011, ApJ, 738, 59 

Rogers, L. A., & Seager, S. 2010, ApJ, 712, 974 

Ribas, I., Guinan, E. F., Giidel, M., & Audard, 
M. 2005, ApJ, 622, 680 

Saumon, D., Chabrier, G., & van Horn, H. M. 
1995, ApJS, 99, 713 

Semenov, D., Henning, T., Helling, C, Ilgner, M., 
& Sedlmayr, E. 2003, A&A, 410, 611 

Seager, S., Kuchner, M., Hier-Majumder, C. A., 
& Militzer, B. 2007, ApJ, 669, 1279 



9 



Stevenson, D. J. 1982, Planet. Space Sci., 30, 755 

Tanigawa, T., & Ikoma, M. 2007, ApJ, 667, 557 

Tcrqucm, C, & Papaloizou, J. C. B. 2007, ApJ, 
654, 1110 

Valencia, D., Sasselov, D. D., & O'Connell, R. J. 
2007, ApJ, 665, 1413 

Valencia, D., Ikoma, M., Guillot, T., & Nettel- 
mann, N. 2010, A&A, 516, A20 

Ward, W. R. 1986, Icarus, 67, 164 

Watson, A. J., Donahue, T. M., & Walker, J. C. G. 
1981, Icarus, 48, 150 

Wuchterl, G. 1993, Icarus, 106, 323 



This 2-column preprint was prepared with the AAS lATJiiX 
macros v5.2. 



10 



